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We present local invariants of multi-partite pure or mixed states, which can be easily 
calculated and have a straight-forward physical meaning. As an application, we derive 
a new entanglement criterion for arbitrary mixed states of n parties. The new criterion 
is weaker than the partial transposition criterion but offers advantages for the study of 
multipartite systems. A straightforward generalization of these invariants allows for the 
construction of a complete set of observable polynomial invariants. 
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1 Introduction 

Quantum mechanical states have a complex description in terms of their density matrix, which 
comprises all information available about a system under a given experimental situation. 
Different density matrices correspond to different states of a system, and allow for different 
predictions on its future behavior. For many purposes, however, we are only interested in 
properties of the state (such as its entropy or purity) which are invariant under unitary 
transformations. 

For systems composed of several parts, or subsystems, there exists a natural tensor product 
structure underlying the state space. For such composite systems, the superposition principle 
gives rise to the phenomenon of entanglement which manifests itself in peculiar "quantum" 
correlations between results of measurements on its different parts 012 G] . To capture the 
essential features of this entanglement, we look for properties of the density matrix that are 
invariant under local unitary transformations, corresponding to a local change of basis in the 
Hilbcrt spaces of the individual subsystems. Such local invariants have attracted the attention 
of people working on the foundation of quantum mechanics and, more recently, in quantum 
information theory 01 |SJ EJ |S] , where entanglement is perceived as a resource for tasks in 
quantum communication and computation. 

In this paper, we present a family of local invariants of a multi-partite quantum system. 
These invariants are derived from an invariant decomposition of the state space of the system, 
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regarded as a real vector space of hcrmitian operators with a scalar product. They have a 
natural geometric interpretation in terms of the length of projections of vectors onto invariant 
subspaces. Such subspaces contain all information either about one local subsystem or about 
correlations between a given set of subsystems. Beyond their straight-forward geometric 
interpretation, these invariants have a number of merits. They can easily be calculated - even 
analytically - for many states, and they are directly connected to measurement data [HI 

The representation of the density matrix as an clement in the real (metric) vector space of 
hcrmitian matrices is well known, and a number of researchers have used a similar approach 
before [111 1121 ^] . Nevertheless, our results go beyond existing work in at least two 
respects. First, the explicit decomposition of the state space into a direct sum of invariant 
sub-spaces makes the identification of invariants quite transparent — it allowed us in fact to 
find a family of new invariants. Second, from the convexity of the set of separable states, 
we are able to derive constraints on the invariants of separable states and propose a new 
entanglement criterion. Here, we will discuss some of its strengths and limitations and apply 
it to a wide class of multi-qubit states. 

Moreover, for each state the obtained invariants are homogeneous polynomials of degree 
2 in the coefficients of a basis decompositions of the density operator into hermitian opera- 
tors. We show how considering polynomials of higher degree one can extend the set of local 
invariants and make it complete UJ. 

2 State tomography 

It is a well known fact that the four Pauli spin matrices <jq = 1,01 — o~xi°~2 = a v ,as = o~ z 
form a real basis of the vector space of the hermitian operators which act on one qubit. 
With respect to the scalar product (A,B) = tr(AB), the basis vectors are orthogonal. More 
generally, for a d-dimensional quantum system, there exists a set of d 2 — 1 traceless hermitian 
generators of the SU(d), which we call ci, . . . cr^-i- One specific choice of these generators 
is the so-called Cartan-Weyl-construction (see, e.g. 0]). Combined with the unit operator 
1 = (Jo, they form a real non-normalized orthogonal basis of the vector space of hermitian 
operators in d dimensions, 

(cri,(Tj) = tr(c7j(7j) = 5 itj d (1) 

Let P = {1, 2, . . . , 7i } be a set of parties and V the vector space of hermitian operators 
acting on the n-partite Hilbert space Ti^ <g> Ti^ ® • • • ® TiS n \ where HS a ' is a Hilbert space 
of the (finite) dimension d a . Clearly, the tensor products of the basis operators form a basis 

B = {4 1] 4 2) • • • 1° < < d 2 a - 1 for all a E P} (2) 

of V. 

Any n-partite density operator p £ V can thus be expanded in the product basis 

1 (!) ( 2 ) («) 

P=~ d 2^ Ciii^-in^h °i 2 •••°i n ' ( 3 ) 

where d = 11"= i d a , and 

^...^tr^^^-'-^O-^^'--^)/ ^ 
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In other words, the expansion coefficients Ci 1 j 2 ...j„ are expectation values of products of her- 
mitian operators. Since these expectation values can, in principle, be measured by local 
measurements (given a sufficiently large ensemble of copies of p) , one can use this method in 
order to determine an unknown n-partitc quantum state with the help of local measurements 
and classical communication (quantum state tomography). 

3 Invariant decomposition of the state space 

Let a = a^a^ ■ ■ ■ er|™' ) be an arbitrary element of the product basis B, and S = {a\i a ^ 0} 
the set of parties, where a acts non-trivially. Using this definition, we call a a S'-correlation 
operator, and the set of all ^-correlation operators Bs- It is clear that B can be written as 
the union of the (disjoint) sets of S'-correlation operators, i. e. 

B=\JB S . (5) 

SCP 

Example: In the case of three qubits, we have eight such sets (with a, b G {1,2, 3}): £>{} ~ {1}, 

#M = K (Q) N = 1,2,3}, B {a . b} = {oj^of \i,j = 1,2,3}, and £ {1 , 2 . 3} = {<r? ) v? ) vj? ) \i,j,k = 
1,2,3}. 

Theorem 1 For each set S of parties, the vector space Vs = span(£>s) is invariant under 
local unitary transformations, which act as isometrics on Vs- 

Proof. Let be a unitary operation which acts on party a. If a £ S, all elements of Bs 
remain unchanged under the action of U^ a \ If, on the other hand, a £ S, then the orthogonal 
set of traceless generators a\ a ^ (i > 0) is transformed into a different set of orthogonal traceless 
generators, i. e. for 1 < i < d\ — 1, 

CT W^) = £(0(C/W)) <£> 

fe 1 

with an orthogonal matrix 0(U^ a ') E SO(d 2 a — 1) ^20]. Obviously, both sets of generators 
span the same set of all traceless operators. □. 

Given a density operator p and a set S of parties, the projection of p onto the subspace 
Vs is given by 

For simplicity, we will often denote the above projection by where the dependence on the 
state p in question will be understood. 

The reduced density operator ps corresponding to the set of parties S can be read off, 

Ps = ps® j-- = E ( 7 ) 

5 S'CS 

where S = P\S denotes the complement of set S. 

By direct substitution into the previous equation one readily checks that the following 
relation between the projection of p onto Vs and its reductions holds 

&G0 = E (-l) |SHS ''/5s'. (8) 

S'CS 
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Due to Theorem 1, local unitary operations rotate a projection £g only within the subspacc 
Vs- Ignoring the normalization constant 1/d leads us to 

Corollary 1 For each set S of parties, the squared length of the projection of a state p onto 
the span Vs of Bs, 

Ls(p)= E (v) 2 p = dtiiis( P ) 2 ) (9) 

is invariant under local unitary transformations. We call L$(p) the S-correlation strength of 
P- 

From Eq.JSJ) we find that S-correlation strength depends solely on the purities of the 
reductions of ps, 

L s (p) = dtii^ s )=dtY(Zsp) = dY / (-V lSl ~ lS ' l te(PS'®^:P) 



dg 



E (-l^H^tr^tr^)] = E (-l) |SHS Vtr04,) , (10) 



da 

S'CS b S'CS 

where ds> = Yl ae s' ^ a an< ^ we nave made use of tr(£g£g/) = dsSs : s'- We notice that the 
S-correlation depends only on the reduced density matrix of the set of parties S, Ls(p) = 
Ls(ps)- Thus, the only invariant which contains information about the total state is Lp(p). 
For pure product states (i.e. tr(p|) = 1 VS) we find 

S'CS aGS' aeS 

which further simplifies to L s nre = 1 for n-qubit systems. 

The idea of local-invariant spaces and correlation strength can be easily extended to a 
scenario where parties can form coalitions, giving rise to different partitions of the set of 
parties. We allow the parties a\ ... a* in a coalition to apply joint operations, or equivalently, 
we treat them as a single super-party acting on higher-dimensional quantum system b. One 
can obtain the required traceless generators for the new party 6 as products of the generators 
of the old parties a\ . . . a^, 

-f, fc =^M: 2) --^r\ (12) 

with (n,...,i k ) ^ (0,...,0). 

Any partitioning can be realized by iteratively joining parties pairwise, say 01,02 —>■ 
b = {01,02}. Using Eq. ljT2|) . one can easily verify that the correlation strength for a set 
S = {b} U S" = {6} U {a M , a Ul . . .} of parties, is given by 

L{b}US' = L{ ai }US' + £{a 2 }US' + £{ ai ,a 2 }US') (13) 

which means that the correlation strengths for coarse partitions are functions of the correla- 
tions strengths of the finest partition. 

A special partition is obtained if we allow all parties to operate jointly as a super-party 
b = P = {ai, . . . , a n }. L{ b } is then invariant under all unitary operations, and thus describes 
a global property of the state. Indeed, we have 

L {P}(P) = E < ct >p - = dtr (P 2 ) - ( 14 ) 
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so that L{py is a measure for the purity of the state p. 

Using Eq. and 03 we can re-write the left-hand side of Eq. El as the sum of all S- 
correlation strengths, 

]T L s =dtT(p 2 ) -1, (15) 

which allows us to state 

Corollary 2 For any state p, the sum of all correlation strengths is given by the purity of 
the state. This implies, in particular, that for states with the same purity, there is a tirade- off 
between local and the different non-local correlations. 

For a pure state p = \ip)(4>\, we have ti(p 2 ) = tr(p) = 1, so that Corollary [2] can be regarded 
as a quantitative expression of the folklore saying that in entangled states, the information 
about the state is contained in its correlations rather than its local properties. 

It is a useful fact that the convex structure of the space V of states is obeyed by the 
subspaces Vs separately, in the following sense: If a state is given by a convex sum of states 
pi, i.e. p = YliPiPi with > for all I and YliPi = 1, then the projection of p onto each of the 
subspaces Vs is the convex sum of the projections of the states pi onto Vs- If p is a separable 
state, it can be written as a convex sum of pure product states. In this case, the projection of 
p onto each of the subspaces Vs is a convex sum of vectors with the squared length L^ UIC , so 
that the squared length Ls{p) cannot exceed L v ™ c . This allows us to formulate the following 
entanglement criterion: 

Theorem 2 Given a multi-partite state p, if there exists a subset S of parties such that 
the S-correlation strength is greater than L^ ure , then p is entangled. 

Corollary 2 then implies that all pure multi-partite entanglement will be detected by this 
criterion. 

It is interesting to note that the strongest entanglement criterion is obtained for the 
finest partition, in the following sense: Let b = {01,02} be a coarsening as in Eq. (|13f) . and 
Ls{p) < L s me for all S C {61, 62} U S' C P. Using (fT^) for the state p and for product states, 
we find 



L{b}uS' — L{ ai }us< + £{a 2 }US' + i{ffli,a 2 }US' 

- rpurc , rpure . T pure . . 

- n { ai }uS> ~T ^{a 2 }US> + ■^{oi,o 2 }US' UOJ 
_ j pure 

— L {b}uS'- 

This means that we do not detect entanglement in any coarse partition, if we do not detect 
it in the finest partition. 

The correlation-strength Ls defined here can be understood as special case of 'g-purity' 
defined by Barnum et aZ.^HI ( see a ls° C3) as the purity relative to a restricted subset of 
observables g. They use this definition to give a 'generalized notion of entanglement': given a 
subset of observables g (not necessarily related to a sub-system), a pure state is 'unentangled' 
iff the g-purity is maximal. Other than using the standard convex roof extension of such 
criterion — for which no closed form is known — , their criterion is limited to pure states. 

The correlation strength criterion is also stronger than the criterion that results from the 
'global' entanglement measure proposed by Meyer and Wallach[TB]. The latter has in fact 
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been proven ^H] to be equivalent to the criterion proposed in applied to multi-partite 
systems. 

For multi-qubit systems, the generalized Bell-inequalities [201 — criterion for the exis- 
tence of a local- realistic description of a set of correlation measurements on a given state — provides 
a weaker entanglement criterion than the correlation strengths. Indeed, if there exists a set 
of local dichotomic measurements on a state p such that a Bell-inequality is violated, than 
it follows that there is a local coordinate system {x, y} such that i =1 iN > 1 (Eq. 

(15) in H2|). Since this is precisely the sum in Eq.© excluding the terms with ik = 3, one 
arrives at Lp(p) > 1. 

However, as one might have guessed by the fact that our criterion only depends on the 
mixedness of p and its reductions, the criterion fails to detect a variety of mixed entangled 
states. This will be seen in the next section where some important classes of multi-qubit 
states will be studied. In fact, it turns out |21| that the proposed criterion is weaker than the 
well known positive partial transposition (PPT) bipartite separability criteria |22| : i.e. any 
state for which a given correlation strength L${p) > 1, will have a bipartite split {A, A} such 
that the corresponding partial transposition results in a non-positive operator p TA ^ 0. 

4 Correlation strengths for multi-qubit systems 

In this section we will compute the local invariants and check the entanglement criterion for 
different classes of relevant n-qubit states. We will see that for most examples the correlation 
strengths can be computed even for large values of n. 

4-1 Dicke- States 

Dicke states |n, m) arc n-qubit symmetric pure states with m excitations (or qubits in state 
10)), 

\n,m) oc^n i |l,.'! l .,l,0,...0) (17) 

i 

with {Ili} being the group of all permutation matrices. We denote p m = \n 1 m){n,m\. Obvi- 
ously, the states \n 1 m) and \n,n — m) will have the same entanglement properties. 

Dicke states appear naturally in quantum optical and condensed matter systems, and have 
also received attention within the field of quantum information (see for example |23 | .|24 | ). 
Due to the large symmetry of these states, the purity of their reduced density matrices can 
be readily computed in the excitation-basis leading to, 




which depends only on the size 1 5*| of the set of parties. 

For m = 1, which corresponds to W states [21], we obtain a S-correlation strength 
^s(p m=1 ) = 7?( n2 + 8|S1 2 — 4(n + l)!^!)- The maximum value is achieved for the corre- 
lation strength involving all parties Lp(p m=1 ) = 5 — — . Entanglement in such states is know 
to be extremely robust to particle loss, exhibiting entanglement even when all but two par- 
ticles are lost. However, we see that using our criterion, entanglement is only detected for 
sets S with more than \S\ > parties. Similarly, for other m values we obtain maximum 
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values that converge as ^ for increasing n and also exhibit a threshold for the minimal size 
of the sets of parties to detect entanglement. 

4-2 Graph and Graph- diagonal states 

Graph-states [251 GUI form a new and very wide class of multi-partite entangled states that 
have been shown to play an important role in quantum information processing: they occur in 
some error correcting codes j^U, provide a universal resource for quantum computation |27| . 
improve frequency standards [22], and allow for secret sharing - Something particularly 
interesting in the context of studying multi-partite entanglement is that they count with an 
efficient description. Given a graph G = (V, E) — i. e. a set of n vertices V connected by edges 
E that specify the neighborhood relation between vertices — graph-states can be conveniently 
defined through the set of commuting observables {K G }2-i, 

K a G = oi a) II ^ foro=l,...,n (19) 

beN a 

where N a denotes the set of neighboring vertices of vertex a ({b : {a, b} 6 E}). The common 
eigenvectors of this set of observables form a complete orthonormal basis of graph-states 
corresponding to the graph G, 

K i G \1>jt) = (-l)' H \1)it)fori = l,...,n (20) 

where the n-dimensional binary array p labels each of the 2™ graph basis states. 

These basis vectors are related by a local unitary = <j%\i/jq), where denotes the 
action of a z on the parties i specified by /i, = 1. 

The stabilizer Sq of a graph-state is the finite group generated by {K G }" =1 . Every element 
in Sq will be a Pauli operator that can be labeled by a binary array v: ap = Ilo=i(-^G 

Thus, in the Pauli operator basis the representative graph state is given by IV'cfXV'ijI = 
1 Sct^gSg ° Vl anc ^ °th er basis states by (V'/jXV'mI = \ J2a es a (~^)^ v<J v- Thus, we notice 
that the jS-correlation strength of a graph states is given by the number of elements of the 
stabilizer that are S'-correlation operators, i. c. Ls(\G)(G\) = \Sq H Bs\. 

By making use of Eq. and the fact that graph states have equally weighted Schmidt 
coefficients for any bi-partition |25| we can also rewrite the .^-correlation length as, 

L S (\G)(G\) = (-l)l s H s 'l2l s 'l- £S '( G ) (21) 

S'CS 

where the Schmidt measure E (G) in respect to partition (S" , S') is determined by the 
Schmidt rank r, E s ' (G) = log 2 (r) gS]. 

An arbitrary n-qubit state can be depolarized by local operations and classical communi- 
cation to a graph-diagonal PHj state. Specifically, this is done by applying sequentially the 
mixing maps T>i(p) = ^(p + K % G pK l G ) for i = 1, . . . ,n. For such a graph-diagonal state pq 
with weights {pp} the 5-correlation strength is given by 

Ls(PG)= E (E^ 1 )^) =PM P (22) 
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with My = E^son^-l)^^ 

This can be readily computed numerically and in some cases also analytically. Below we 
discuss the class of GHZ-diagonal states. 

GHZ-diagonal states: The generalized GHZ-states correspond to 'star'-graphs |25| . 
The representative graph state is given by \ip^) = -i=(|0) 2 |0 . . . 0) x + \l) z \l . . .l) x ). The 
remaining GHZ-basis states are \ipp), where the first bit of the binary array p, determines 
the phase, while the remaining n — 1 bits correspond to the bit-flips in the n — 1 last qubits. 
Here, we will label each graph-state basis vector either by the the binary array p or by the 
corresponding integer in the decimal basis p £ [0, 2" — 1]. 

Since for GHZ-states the Schmidt measure is E s {\GH Z)) = 1 for all partitions with 
exception of S = {} and S = P for which it vanishes, one can readily show that (Ea l22fl 

L P {GHZ) = 2"- 1 + S n , ovcn and L S {GHZ) = d>|, cvcn for S C P. (23) 

Following the previous notation for the stabilizer group elements, the set of Pauli operators 
cr,7 £ S G r\B P is given by „4 odd = {V} = {(1, x)\x = 0, . . . , 2™" 1 - 1} for n odd, and „4 cvon = 
„4 odd U {(0, 1, . . . , 1)}. The P-correlation strength for GHZ-diagonal states can then be shown 
to be (EqEU, 

2 «-i_! /V- 1 -! \ 2 

L P {p)^2 n - 1 (PM-PM+2-0 2 + ^,even ]T (p M + P M+2 n-x)(-l)^ , (24) 

where we have used ^2 S (— 1)' T = bp$- n - According to our chosen convention, the state 
corresponding to the index-label /j, + 2 n ~ 1 is the same than that with label \x but with opposite 
phase (±). Hence, the weights and of GHZ states with opposite phase always 

appear together as a difference (and also as sum in the even n case). Since the reductions of 
a GHZ are always separable, it is clear that Ls(p) < 1 for S C P. 

These expressions can be further simplified if we restrict to a still very relevant class that 
can be obtained from a general state by a further local depolarization step [32]. The class is 
fully specified by 2™ _1 parameters, 

Pi = P 4 +2™-i for i > 1, and A = p - p 2 «-i. (25) 

Such states play the role of generalized multi-qubit Werner-states and have the nice 
property that the positivity of the partial transposition in respect to a subset A is easily 
checked: p TA > <^> A < 2p^ (where jl is the binary array with pi = 1 iff i £ A), 

For such class of states the P-correlation (Eq. I24J) only depends on the parameter A (for 
odd n). 

L P (A) = 2"- x A 2 + 5„, cvcn ^2 J2 + A j (26) 

In the remaining of this section we use the previous results to compare the performance 
of the correlation strength with other entanglement criteria: 

Lp{p) -> P > (2"" 1 + 5„,even)"5 

For noisy GHZ: p\GHZ)(GHZ\ + (1-J>)^1 { Bell Ineq. -► p > (2"" 1 )"5 

NPPT p > (2 n - 1 + I)- 1 
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• L s (p) > 1 <f> (GHZ\p\GHZ) > \ => N-distillability. That is, the correlation strength 
criterion is neither stronger nor weaker than the MES-overlap criterion. 

• The correlation strength criterion is neither stronger nor weaker than realignment cri- 
terion [331 . It is known that the realignment method can detect some PPT bound 
entangled states, while the correlation strength can not. On the other hand, there are 
simple examples — even in the 2x2 case (see j!2|) — , where the entanglement is not 
detected by the realignment criteria, but £{1.2} > 1. 

( L P {p) -» N>8 

• It detects NPPT multi-partite 'bound' entanglement \6h\ . < Bell Ineq. — > N > 8 

[ NPPT -» N > 4 

• If p has m or more positive partial transposes Lp(p) < Lp(A c ) where A c = m * 1 . 

• k-separabihty: If a GHZ-diagonal state can be written as p = J2iPiPi ® • ■ ■ ® Pi there 
will be m = 2 fe ~ 1 — 1 PPT, and one arrives to Lp(p) < 2 n ~ 2k+1 . This coincides with 
the strongest upper-bound obtained from quadratic Bell-Inequalities |36j in restricted 
experimental settings |37]. Using the property L$(pa ® Pb) = L(pa)L(pb) it is also 
possible to derive k-separability criteria for general states |21j . 

5 Towards a complete set of invariants 

The local invariants L$ do not form a complete set of invariants, i. e. they do not contain 
all information about the entanglement properties of a given state. However, the formalism 
used in this paper allows us to identify a larger class of invariants, many of which also have 
a straight-forward geometrical interpretation. 

From the proof of Theorem^ it follows that the transformation properties of the subspaces 
Vs are closely related. In order to show how this can be used for the construction of invariants, 
we first define the S'-correlation tensor C5, which is composed of the components of the 
projection £g in |jf)J), 




One can easily see that a contraction of two such tensors with respect to a index i v at the same 
position is invariant under local unitary operations, i. e. under orthogonal transformation O 
which affect this index: 




— ^ ] ...^i' v i'l^...i'^ ... (28) 

i' v 

Any complete contraction of correlation tensors, i. e. a polynomial in the expansion coefficients, 
in which indices are either zero or summed up pairwise, is thus a local invariant. Examples for 
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such polynomials are CfyfcCyoCiofc, cojooCjjOiCij' koc C j' ki (where, as usual, the sum is taken over 
all indices which occur twice), the correlation strengths L$, and other objects which can be 
interpreted as scalar products, such as the scalar product of p ai a 2 a 3 with the tensor product 
of P{ai},P{a 2 } and p {Q3} , 

(P{ ai } ® P{a 2 } ® P{a 3 },P{ ai a 2 a 3 }) = ^ ^OOCOjOCoOfcCyfc. (29) 

i,j,k>0 

Unfortunately, it is not possible to construct a complete set of local invariants using the 
construction above: already for the case of two qubits, there are seven independent invariants 
which can be written as contraction of correlation tensors; the two remaining invariants are 
functions of the determinant and sub-determinants of the correlation tensor |13| . 

Following the ideas of we can can generalize the above constructions in order to obtain 
a complete set of polynomial invariants in the coordinates Cjj...i n : 

Any homogeneous polynomial / of degree fc, i.e. 

/(c»i...i„) = E M i\,...,ii,il ** Cil-iiCi?...** '"Cif ...«*> ( 30 ) 

is invariant under local unitary transformation, iff the corresponding observable 

M = E (31) 

i 1 ,...,i n 

on <g) . . . ® TL {n) ) m commutes with all unitaries of the form {U^ ] <g> . . . <g> [/ 1 l™ ) )« ,fc . For 

example, the proposed invariants Ls are homogeneous polynomials of degree 2 for which the 
corresponding observable is Ms = X^o-eBs 17 ® G - ^he ma P ^ s a vector space isomorphism 
from the algebra of observables M on H m , with [M, (C/f 5 <g> ... ® C/l" ) )® fe ] = for all 
local unitaries, onto the the algebra of invariant homogeneous polynomials of degree k in 
the coordinates Cjj...j n . A measurement of the observable M on k copies of the state then 
evaluates the corresponding polynomial invariant M, i.e. fM{ci 1 ...i n ) = (M 7 p® k ) = tvMp® k . 

Note that the orbits of the compact group of local unitaries can be separated by some 
polynomial invariants (see |38| . p. 133). This means that two states p\ and pi can be trans- 
formed into each other by local unitaries if and only if the evaluations of all polynomial 
invariants / coincide. Moreover, the algebra of invariant polynomials is finitely generated and 
the generators can be chosen to be homogeneous polynomials. In order to find a complete 
set of invariants, that generates the whole algebra of polynomial invariants, it is therefore 
sufficient to find the generators of the algebra of invariant homogeneous polynomials for dif- 
ferent degrees k, or equivalently the algebra of hermitian matrices M, that commute with 
{U^ ® . . . ® U n n ^)® k for all local unitaries <g> . . . <g) Un 1 ■ The proposed set of invariants 
{Ls} arc then an example of generating set of homogeneous polynomials of degree k = 2. 

In [H] the construction of a set of (possibly non-hermitean) generators F can be found. 
From it, a set of hermitean generators can be obtained by symmetrization: M\ = F \ F and 
Mi = F ~2^ ■ For each generator M the coefficients M^ ^ ^i q ^k then can be computed 
from Eq. I|31|l . which determines a set of generators /m for the invariant ring of degree k by 
Eq. JU ■ 
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6 Conclusions 

In this paper, we have investigated local invariants in terms of correlation operators. The 
related geometric interpretation of the space of coherence vectors and its LU-invariant sub- 
spaces allowed us to explicitly write down a subclass of such invariants. Using a convexity 
argument, we have given a simple entanglement criterion for multi-partite states. 

We have shown that the criterion fails to detect important classes of entangled states 
— including some NPPT states. Nevertheless, we think that our criterion is of interest. The 
criterion is simple and with a clear interpretation. It is easily computable even for large num- 
ber of qubits and some closed forms can be obtained. In particular they are more efficiently 
computable than the PPT-criteron that involves a matrix diagonalization. This might be 
especially usefull in the study of entanglement in many-body systems. On the other hand, 
Lp gives information on all partitions so that one does not need to check the positivity of all 
2™~ 1 — 1 possible partial transpositions. 

From the experimental point of view our criterion also offers some advantages. Firstly, 
the proposed criterion are state-independent: it does not require previous knowledge of the 
state nor optimization over possible experimental-settings, like in the case of Bell-inequalities 
or general entanglement witnesses. To measure the S'-correlation Ls(p) efficiently one can 
measure a single observable Ms = J2aeBs ° ® a on ^ wo C0 P^ es P ® P requiring only bi-local 
operations between equivalent qubits in each of the copies. Alternatively one can opt for a 
more straightforward approach and measure the correlation functions (o'i^o'i^ • ' ' "i™^ > as 
one does in quantum tomography, with the advantage that entanglement might be detected, 
i.e. Ls becomes larger than one, long before the complete tomography is over. 

Finally, it remains an open question whether the values achieved by the correlation 
strengths are related to the distillability properties of the state. 
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